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Using the scattering transform for n*^ order linear scalar operators, the Poisson 

bracket found by Gel'fand and Dikii, which generalizes the Gardner Poisson bracket for 

2 ' the KdV hierarchy, is computed on the scattering side. Action-angle variables are then 

constructed. Using this, complete integrability is demonstrated in the strong sense. Real 

action-angle variables are constructed in the self-adjoint case. 
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^vj ' 1. The Gel'fand-Dikii hierarchies. 



In this paper we construct action-angle variables and show complete integrabil- 
ity in a strong sense for the Gel'fand-Dikii flows, which generalize the KdV hierarchy. 
Action-angle variables were first constructed for the KdV hierarchy in 1971 by Za- 
kharov and Faddeev [ZF]; and for the nonlinear Schrodinger hierarchy in 1974 by 
Zakharov and Manakov [ZM]. McLaughlin [Mc] derived action-angle variables for 
the Toda lattice and Sine-Gordon equations. 

In a previous paper [BSl] we constructed action-angle variables and showed 
rS ' complete integrability for a general class of Hamiltonian hierarchies based on first 

c^ ■ order n x n isospectral operators. It was shown that in some cases, e.g. the three 

wave interaction, it was necessary to introduce flows nonlinear in the scattering data 
in order to obtain the maximal number of commuting Hamiltonian flows required 
for complete integrability. This is in contrast to the KdV and nonlinear Schrodinger 
equations in which the necessary commuting flows can be taken to be linear in the 
scattering data. Consequently, the method of Lax pairs does not in general yield 
enough commuting flows for complete integrability. In the case of the Gel'fand- 
Dikii flows, however, the Lax pairs themselves do yield enough commuting flows for 
complete integrability, and these flows are linear in the scattering data. 

Gel'fand and Dikii [GD] constructed a hierarchy of Hamiltonian flows which 
generalize the KdV hierarchy based on the Schrodinger operator L = D^ + u. They 
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considered flows based on the n*'^ order scalar differential operator 

" Id 

L = y^Uj{x)D^; D = --—; Un = I, Un-i = 0. 
^-^ I ax 

where Uj = Uj{x), j < n — 1, are elements of the Schwartz class of functions S{R). 
The corresponding flows are given by 

^ = [4/",L], fc = l,2,... (1.1) 

where k/n is not an integer, L^' "^ denotes the formal n*^ root of L considered as a 
pseudodifferential operator with symbol 

Kx,i) = ^Wj{x)C^; 

and B^ denotes the differential part of a pseudodifferential operator S^. In the 
case n = 2, L is the Schrodinger operator, and (1.2) is the well-known Korteweg- 
deVries hierarchy (k odd). For n = 3 and k = 2 equations (1.1) yield the Boussinesq 
equation. 

Gardner [G] showed that the KdV equation is Hamiltonian with respect to 
a natural Poisson bracket. Gel'fand and Dikii showed that the flows (1.1) are 
Hamiltonian with respect to the Poisson bracket 

r,s 

for the n^^ order case, where i is the matrix of differential operators 

n-l-r-s , , 

^rs= Yl [{ )nr+s+k+lD'^ - { ){-D)'^Ur+s+k+l] (1-3) 

fc=0 

(We have introduced the factor i in (1.2), by comparison with [GD], so that the 
bracket is real in the self-adjoint case L = L* .) We assume throughout that the 
Ujjj < n — 1, belong to the Schwartz class 5(i?), so the terms {—D)^Uk+r+s+i can 
be integrated by parts. After some manipulation, (1.2) can be written 

{F, G} = 

J -r^,^-^ ^-^ ^-^ '^ S 6u^-r-l SUs SUs 5u^-r-l 

^ °o 7=1 r=0 s=0 7T-1 s s 7ri 



GELFAND-DIKII FLOWS 3 

We will need this form below. 

The scattering transformation of (1.1) as developed in [BDT] is based on the 
wave functions -i/; which satisfy Li/; ~ z'^if). This v}^ order scalar equation can 
be written as a first order system in the standard way as follows. Let Cjk denote 
the matrix with a 1 in the jk place and zeroes elsewhere. Let V'l be a solution of 
the scalar equation; the column vector -(/' = ('*/'i7'*/'2 • • •V'n)* satisfies the first order 
system [D — J^ ~ q)i^ = 0, where 
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We denote by P the linear space of all such matrices q with entries in S{R). Let $ 
and \1/ be n X n matrices which satisfy the first order systems 



D^ = {J, + q)^ 



(1.5a) 



L>^ = -^(J^ + g) 



[1.5b) 



These wave functions are uniquely determined by specifying their normalization 
as a; ^ — oo; that normalization is given in (1.9) below. Define the matrix a by 
a = \l/$; it is easily seen from (1.5) that a is independent of x. 

The entries of the first row of $ satisfy the equations L(p = z'^cf) while those of 
the last column of \E' satisfy the equation L^'^ — z"^"^, where L^ is the transpose of 
the operator L: 



L^^ 






-D 



)^Wj* 



Moreover, (1.5a) shows that $ is the Wronskian of the functions of its first row. 
For n > 3 let E be the set of rays in the complex plane given by 

T, = {z E C : Re ia-' z = Re ia z for some < j < k < n} 



where a = exp27ii/n is the primitive n^^ root of unity. We label the rays Ej, j = 
0, 1, . . . 2?i — 1, with Eq the negative imaginary axis, and the numbering proceeding 
in the counterclockwise direction. We label domains of the complement C\E by 
Oo7^i7 • • • • (When n = 2 the appropriate set is E = i? and minor bookkeeping 
changes must be made in the following analysis; since the results were known for 
n = 2 we omit the details.) 
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Fig. 1. Fundamental Domain 

It will sometimes be convenient, as in [BDT], to work in a local ordering, that is, 
in an ordering that depends on z. We define the j^^ ordering of the roots to be 
{«!, a2 . ■ ■ ctn}, where 

Re iaiz > Re ia2Z > ■ ■ ■ > Re ianZ for 2 G Oj. (1-6) 

The ordering in Oq is given by {1, ct, ct~^, a^ . . . } while that in Oi is given by 
{1, a~^j ct, a~^, . . . }. The rays Eq and Si are given by 

Eo = {^ /m e < 0, Re ta^C = Re ia-^i, J = 1, • • • , [-]}, 

ft 

T.^ = {i:Imi<Q, Rei a-^i = Re i a^+'C. J = 0, ...,[-] - 1}. 

The transformation of a matrix from the j^^ to the (j' + l)^* ordering is accomplished 
by conjugation by a permutation matrix ttj. We see that tvj = 7rj_|_2, and that 
the structure of the permutations is different for even and odd n. The associated 
permutations are 

TTo = (23) . . . (n — 1, n) and tti = (12) ... (n — 2, n — 1), n odd; 

and 

TTo = (23) . . . (n — 2, n — 1) and tti = (12) . . . (n — l,n), n even. 
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We define the diagonal matrix J{z) or Jj to be the diagonal matrix with entries 
(cti, q;2, • • • ctn) where the roots are in the j^'^ ordering. Similarly, let 



A, 



1 1 ... 1 

tti a2 • • • ttn 

^1 ^2 • • • '^n 



where the roots are in the j^^ ordering; and for z G Oj, define A^ = d{z)Aj, where 
d{z) = diag{l, z,z'^, . . .z'^~^). Then A j A* = n/ where / is the n x n identity matrix, 
and Jz^z = zAzJ{z) for each z G C\E. 

We now make the gauge transformations $ = A;j(/), \1/ = 'i/jA~^, qAz = AzQz to 
get the equations (J = J{z)) 

D(P = {zJ + qz)(p, (1.7a) 

D-^ = --^{zJ + qz). (1.76) 

Recall that a{z) = \l/$. It is clear that a{z) is invariant under the gauge trans- 
formation 9? ^ $, f/j -^ \1/; hence a{z) — 'i(^{x,z)4>{x, z). In [BDT] fundamental 
solutions of systems (1.7a,b) are constructed in each component of C\E by a tech- 
nique involving Volterra integral equations for successive exterior products of their 
columns. The solutions are constructed to satisfy the asymptotic conditions 

lim (Pe''''^ = l, limsup||(/)e-^^^'^|| < +oo; (1.8a) 

lim e'^^^^ = 1, limsup ||e'^^^*|| < +oo, (1.86) 

where J = J{z). The fundamental solutions $ and \E' satisfy the asymptotic condi- 
tions 

lim Aj^$e-^^^^(^) = 1, lim e^^^-^^^^^A^ = 1. (1.9) 

X — » — oo X — yoo 

As a; ^ — oo, respectively -|-oo, $ and \1/ are asymptotic to U = A^e*^^"'*-^-' and U~^ 
respectively; a simple calculation shows that 

Ujk = L>^-le^^"'=^ {U-^)jk = (-L>)'^-'=-^e-^^^"^. (1.10) 

Let E j be a ray of E and let ^ belong to Ej . Denote by (p± the boundary values 
of the wave function (p defined on the components of C\E separated by Ej. In 
[BDT] the scattering data for L associated with Ej is defined by the jump relations 
"^(0 = 4'- (p+'^j- III general, the scattering data consists of discrete data, as well, 
which lives on those points z for which the construction of (p and t/j breaks down. 
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We assume in this article that there is no discrete component of the scattering data, 
so that L is uniquely determined by its continuous scattering data. 

As in [BC] and [BSl] we may instead define a scattering matrix for the first 
order system (1.7a) as follows. Fundamental solution matrices which satisfy the 
asymptotic conditions (1.8) and (1.9) can also be constructed for ^ G E\0. For 
^ G E, let n^ denote the projection 



{'n^s)jk 



Sjk if Re iC{Xj — Xk) = 
otherwise. 



Then the limit 



s(0 = lim n^e-^^^^(/)(x, = lim U^e''''^-^ A7^<5{x,0 (1-12) 

X — > + oo X — > + oo ^ 

exists. The matrices s(^), as ^ ranges over the rays Ej, determine the continuous 
data V, and, in the absence of discrete data, v determines s. In §3 action-angle 
variables for the Gel'fand-Dikii flows are constructed directly from the entries of s, 
not V, so it is important to pass to s. 

Lemma 1.1. s = a = tpcj). 

Proof: Writing (j) = me^^^'^ and i/j = e~'^^^'^w, the functions m and w are bounded 
on — oo < a; < oo by (1.8) and (1.9). Therefore wm — e'^^'^a e"*^^"^ is bounded for 
X E R, and so Il^a = a. Now 

s(0 = lim U.e-'''^-^ (p = lim n.e-^^^-^-^'^o 

X — »00 X — >CXD 

= lim U^e-'^'^-^W-^e'^'^-^a = lim [U^e-'''^-^w-^e'''^'^]a 

X — >CX) X — »oo 

since 11^ (Sa) = {Il^B)a for any matrix B whenever a is in the range of H^. As 
X tends to oo, w~^ tends I; and since H^ projects onto the oscillatory part of 
exp{ix$,adJ{^)}w, 

lim n^e-'^^-^w-^e^^^^ = /. 

X — >oo 

Hence s = a and Lemma 1.1 is proved. 

Since a lies in the range of the projection H^ its only nonzero entries occur 
where Re iaj^ = Re iakC- ^^ the local ordering, a(^) has a block diagonal structure 
consisting of diagonal entries or 2 x 2 diagonal blocks. We label the 2x2 blocks by 
B^, p — 1, . . .n — 1 as indicated. The block diagonal structure of a(^) is different 
for even and odd n as indicated in the following diagrams: 
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/an 



B2 




B4 



n odd 



Bn-l / 



S, 



\ 



(^nn ' 



Si 



/an 



So 



S, 



n even 



n-2 



/ 5l 



s 



n-1 



Si 



s„ 



a,. 



tti^+l,;y 



(i'v,v^\ 



O-v+l^v+l 



detS^ 



Fig. 2 Block Diagonal Structure 

Lemma 1.2. The wave functions $ anc? \1/ and scattering data a(^) are invariant 
under the rotations z -^ az; in particular, a(^) = a(a^). Consequently the scat- 
tering data is determined entirely by the matrices a(^) for $, on the rays Eq and 

El. 

In Lemma 1.2 it is understood that we are working in the local representation. 
The invariance of the wave functions and scattering data under rotation by a means 
that all data is contained in the wave functions defined on the fundamental domain 
in Fig. 1. 

Proof: First note that J^, A^^, and zJ{z) are invariant under the rotation z -^ az, 
since the matrix representations are local. Consequently, equations (1.5a,b) and the 
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matrix U of (1.10) are invariant. By uniqueness, the wave functions $ and \1/ are 
invariant, and therefore so too is a(^) = \l/(a;, ^)$(x,^). 

Later we show that the scattering data possesses an additional symmetry when 
L is self-adjoint. 

2. The Poisson brackets. 

We now turn to the computation of the Poisson brackets of the scattering 
variables. The entries of the scattering matrix a(^) are regarded as functionals on 
P, and we compute their corresponding Poisson brackets. As in [BSl], the gradients 
dttjkiO/^Q do not decay as |a:| -^ oo. Therefore the integral in (1.2) does not have an 
absolutely convergent integrand, and it is necessary instead to use the regularization 

hm ^ r y[lJ-^]^-^lp^dx. 

TV^oo _AT ^^-^ OUs OUr 

r,s 

This limit exists only in the sense of distributions; so the precise interpretation of 
the calculation is as follows. Given any pair of functions /, (7 G T>(Tj), define the 
functionals 

F{q) = / ajkiOfi^dC, G{q) = I aim{ri)g{vi)dr]. 

JT, JT, 

Then the formal calculation 

Y^[i ] dx 

N^^-^ OUs OUr 

r,s 

'^^ Jt.Jt.J-N~' OUs OUr 

= {ajk{i),ai^{r])}f{i)g{rii)didr] 

JT JT 

furnishes the defining relation for the distribution {ajfc(^), 0^^(77)} G V'iTi x E). 

Let (5$ and 5a denote variations in the wave function and scattering data due 
to a variation in the potential matrix q (i.e. in the coefficients of L). The following 
identity may be established by replacing qhj q + eSq, forming the equation for the 
difference in the wave functions, and taking the limit as e — * (cf. §D of [BC]): 

D(*5$) = ^5q$. (2.1) 

Integrating this identity over (—00, 00) and using the fact that 5(j) and St/j vanish at 
—00 and +00 respectively, we obtain 

/.oo 

da = i "^dq^dx. 



N- 
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Therefore 



tr[^6q^ekj]dx = ^ / tr[^ekj'^en,i+i]Suidx. 

-CX) 7 Q J —OO 

Hence 

(-0%^ = ($efc,*)z+i,„ = $z+i,fc(x,e)^,-„(a:,0 = (I^Vfc)V'j (2-2) 

where (/)a; and tj^j are the entries of the first row of $ and last column of \1/. The 
Poisson bracket between entries of the scattering matrix a is therefore (using (1.4)) 



oo n 7-1 r 

-oo -y=l r=0 s=0 ^^•'^'' 



Here -i/'j and (/)fc are evaluated at (x, ^) and -i/;; and (jij^ are evaluated at (x, ry). 

When ^ and tj are on different rays we must choose a common representation 
for the matrices a(^) and a{ri); in the following, we work in a fixed representation. 

The integral in (2.3) is evaluated by expressing the integrand as an exact de- 
rivative. More precisely, 

Lemma 2.1 Let w-y,7 < n — 1, belong to the Schwartz class. Then 

-i ("^ 
{a,fc(0,az^(^)} = T;rr^ / D[W{ct>k{O.Uv))W{ct>m{v),^M))]dx, (2.4) 

S 'I J — oo 

where 

n 7— 1 
7=1 r=0 

Proof: For scalar functions / and g, 

gLf-fL^g = DWif,g). 
In particular, if Lf = ^'^f and L"^ g = rj'^g, then 

DW{f,g) = {C-vnf9- (2.6) 

It follows from (2.5) and (2.6) that 
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= (j)ki^lW{(j)m, Ipj) - (j)mtpjW{(j)k,1pl) 
n 7 — 1 

7=1 r=0 

Since we assume the m^, 7 < n — 1 , are in S{R) this expression is equivalent modulo 
an exact derivative to 

n 7— 1 

7=1 r=0 

Comparing this with (2.3) we get (2.4). 

We now evaluate the integral in (2.4). As noted in the beginning of this section 
the integral in (2.4) must be interpreted in the sense of distributions. We obtain 



w{MO.Mv))w{c^m{v).i^jiO) 



{a,ki0.aimiv)} = - lim " -.v.., ..v...-- v...w., -..^ 



N^oo ^"^ - 7]"^ 



N 



-N 



where (t>k{0 denotes (/)fc(x,^), etc. Since Un = I and the remaining coefficients Uj 
tend to zero rapidly as a; ^ ±00, 



n-l 

|a;| ^ 00. 

s=0 



Vr((^»~ Woi(P,ij) = J2iD^-'-^(P){-DyiJ as 



From \l/$ = a and (1.10) we have, as x ^ —00 and x -^ +00 respectively, 



and 



«r- .-^C^, , tt 



^.-E«-(0-7;:^e-^- t 



n^'^-i ■" n^ 



■? ^fn-l 



:Z -,-«2;^"i 



From the block diagonal nature of a(^), the sums for (pj and t/jj extend only over 
the one or two element set {r : Im ^(a^ — ctj) = 0}. 
Straightforward calculations show that 

en _ ^n 

and that as x — ;► ±00 

^-L_w{MO. Mv))w{<Pm{v). MO) 
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is asymptotic to 






and 



Pi' 

respectively. Hence 

{ajfc(^),az^(?7)} = -^apk{Oarni{v)Rpi,rji^^v) - aip{v)ajr{ORkp,mM^v) (2-9) 

where 

( — D^i^^ — ri^) „ix(£,{ap-aj)+r]{ar--ai) 

R^ {^^r]) = - ^7^:^ — T^aia^ lim — — — . (2.10) 

Lemma 2.2. For ^, ?7 in Eq U Ei ^/le limits in (2.10) are 






(_l)n-l 

± m5ip5jr — TT^T^'^^ ~~ V)sgn{l - j) 



(2.11) 



Here d{$, — rj) denotes the complex valued distribution defined with respect to the 
complex measure d^ on the rays of E such that 

f{vnC-v)dv = f{0- 

E 

Brackets for general ^, r] are obtained using the symmetry in Lemma 1.2. 

Proof: We use the O*'* representation on both rays. Note that the block diagonal 
structure of a on Ei is the same in either the O*'^ or the 1** representation and so 
is the same as that shown in Fig. 1.2. 

The first term in (2.11) arises when ^{ap — aj)+r]{ar — cti) vanishes identically in 
^ and ?7, that is when j = p and r = 1. In that case the argument of the exponential 
vanishes and we get a simple principal value term. Note that 

tn _ „n 

^ ^ (2.12) 

^ttj - Tjai 
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is regular, since aj and ai are roots of unity. 

Due to the block structure of a(^), the sums in (2.8) extend only over the sets 
where Im ^ {ap — aj) — Im ry (ct^ — ai) = 0; hence the exponential term in (2.10) 
is purely oscillatory. We use the identity 

lim p.v.^ = ±ni sgn{X)6{C-r]), X e R\0, (^ - r]) e R, (2.13) 



X — >±oo 



to evaluate the limit when the argument of the exponential does not vanish. Since 
(2.12) is regular, the limit in (2.10) vanishes by the Riemann-Lebesgue lemma when- 
ever only one of the factors {$,ap — rjai), {^aj — rjar) vanishes. Since Sq and Ei are 
on different orbits under multiplication by roots of unity, these factors can vanish 
only if ^ and r/ are on the same ray. If both factors vanish, then j — r and p — I 
and we must evaluate 

^ix(S,-r])(oii-Oij) 

lim 

when Im^^ — r]){ai — aj) = 0. In evaluating this limit, we must keep in mind that 
it is defined in the sense of distributions; so we really need to consider the limits 

r gia;(C-'?)(«j-«i) 



X — >±oo 



X 



r Qix{S,-ri){aj-oti) 

lim / 77 77 sfiv)d{aj - ai)r] 



where / is a smooth test function on the ray E^., r = 0, 1. Making the change of 
variables s = ^{aj — ai), t — r]{aj — ai) we get 

psgn{l-j)oo ix{s-t) f 

- lim / /( )dt 

x^±cx) Jo ^-t OLj - Oil 

= ±7ri/( )sgn{l - j) 

aj — ai 

= ^T^ifiO^gnil-j). 
Therefore in the sense of distributions the above limit is 

gix(^-T;)(a;-aj) 

lim =±ni5{^-7])sgn{l-j). 

x^±oo t fj 



Equation (2.11) follows, since 



en _ r 

^ 'I SO- ..\ ../-n-1 



C-v 
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Using lemma 2.2 and (2.10) we get 



(-1)" 

+ (-l)'"7rza;fc(0ajm(0 cn-i ^^^ ~ v)[s9n{l - j) + sgn{k - m)] 
where 



(2.14) 



Caj - rjai ^aj - rjai 

In (2.14) all matrices are in a fixed (global) ordering. We shall always take the 
ordering to be that of Qq. For n = 2 (2.14) gives the known result for the KdV 
case. (cf. [NMPZ] p. 41; the factor of 1/2 is due to different normalizations of the 
Gardner Poisson brackets.) 

3. Action-angle variables; general case. 

The construction of action-angle variables for general first order nxn isospectral 
operators was given in [BSl]. While the situation for the n^'^ order scalar case is 
similar, there are some differences. The local part of the bracket (2.14), i.e. the 
part containing the delta term, is the same, apart from the factor {n^^~^)~^ , as in 
the case of first order systems. The nonlocal term, - involving p.v.{^aj — r]ai)~^ 
- is more complicated than that in the systems case. Finally, because of the block 
structure of the scattering data, we need only consider the case of 2 x 2 matrices, 
for which the structure of the action-angle variables is especially simple. 

It is convenient to use the positive real, rotationally invariant measure C~^d^ 
on E. We denote by S{^/r]) the associated distribution defined by 

E ? 

A simple computation shows that S{^/rj) = ^ 5(^ — ry), the latter being the distri- 
bution defined in the previous section. 

We label the 2x2 blocks of the scattering matrix a hj u = 1, . . .n — 1, as in 
Fig. 2, and define as the canonical variables associated with the u^^ block, 

PA0=^{-0 log , g^^ =—- log -— , 3.1 

Aj, 47rz a^^^+i{^) 

where A,^(^) = ai,i,a^+i^^+i — ai,^i,+ia^+i^i, is the determinant of the 2x2 block Bi,. 
In (3.1), ^ is on Sq or Ei according as u is even or odd; and we work in the local 
representation. The phase factor a'^ inserted in the definition of the Qi, is needed 
to get the right symmetries in the self-adjoint case, discussed in the next section. 
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It simply involves changing q^, by an additive constant, and so does not affect the 
Poisson brackets. 

Theorem 3.1. The variables p^.q^, satisfy the canonical commutation relations 

with respect to the positive measure ^~^d^. 

Proof: The distribution-valued Poisson bracket (2.15) is a sum of two terms, a 
local term, involving the distribution d{^ — r]), and a nonlocal term, involving the 
principal value terms such as p. v. (^ctj —r}ai)~^. (Note however, that (^ctj — r]ai)~^ 
is singular only if ^ and rj are on the same ray, since Eq and Si are on different 
orbits under multiplication by a root of unity.) As in [BSl] we denote the local part 
of the Poisson bracket by [ajfc, aim] and the nonlocal term by < ajk, aim >■ 

Recall that the relations (2.14) refer all quantities to the O*'^ representation, 
whereas the definition of the canonical variables in (3.1) uses local representations 
on each of the rays. To convert the representation from 0^^ to 1**, we conjugate 
by the permutation ttq; and this interchanges a,^^, with a,^_|_i^,^_|_i and a^^^^i with 
ai^j^iy, so Pi, is unchanged and qi, is replaced by {const. — qi,). Since the local 
component of the Poisson bracket vanishes when ^ and ry are on different rays, none 
of the following computations are materially affected. 

By (2.14) we have for all j, k, ^, r] (with j < k, say) 

(-1)" 
{ajjiO.akkiv)} = ajjiOakkiv) ^2/£ ^n-i iQjM^^ v) - QjM^^ v)] 

1 

+ (-l)'"7riajfc(0afcj(0^^m'^(^ - v)[s9n{k - j) + sgn{j - k)] = 0. 

nq 

The first set of relations in (3.1) follows if we can show that {A^, A;^} = {A^, ajj} = 
for any fx, z/, and j. Note that for the local part of the bracket, [ajfc(^), aimiv)] = 
whenever ^ and r] are on different rays, or whenever ajk is strictly "northwest" or 
"southeast" of a^^n, due to the factor [sgn{j — l) +sgn{m — k)]. The same result can 
also be deduced from the block diagonal structure of a. Thus [A^, A,^] = [A^, ajj] = 
0. Let us consider the nonlocal bracket < A^,ajj >. We have, by the derivation 
property of the Poisson bracket (cf [BSl]), 

< Ai,,ajj >—< a,y;^a,y_|_i^;y_|_i — a,^^,^+ia^-^.i^iy,ajj > 

= tti^jy < a;y-(-i,,y-(-i, ajj > + ttjy-i-i^jy-i-i < ajyjy, ajj > — aj^^jy+i < aj^^i^^, ctjj > 



v?{^r]) 



■:—[(^u,u+lO'u+l,uO'jj[Qu+l,j ~ Qu,j + Qu,j ~ Qu+l,j] — 0. 
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A similar calculation shows that < A^, Aj^ >= 0; and the first set of commutation 
relations in (3.1) is proved. 

Next, 

^log-HdhH^log^^+h^j = {loga^+i,^ - loga^,^+i, loga^+i,^ -loga^,^+i} 

1 . 1 

We claim that the local part of the above Poisson bracket vanishes. It certainly 
vanishes if ^ and ry are on different rays. If ^ and ry are on the same ray and ^ < u, 
the local component vanishes due to the block structure of a; while if n = u then 
the first and fourth terms vanish while the second and third cancel. 

Therefore, the above Poisson bracket reduces to the nonlocal terms, viz 

(-1)" 



— {Qn,u+i — Q^+i,v) + Qfi,u — Qfi+i,u+i] = 

and {quiO^Q^iiv)} = 0. 

Turning now to the computation of {pu{0^ ?m(^)}' ^^^^ note that the local part 
vanishes when n ^ u. As before, this follows from the block diagonal structure of 
a or from the structure of the factor [sgn{j — /) + sgn{m — k)]. When n = u and ^ 
and r] lie on the same ray, 

[loga^^(0,loga^+i,^(?7)] = [a^^, a^^+i,^] 

We leave it to the reader to verify that the other three terms give the same result, 
and that [A^(^), g^(?7)] = 0; hence 

[pAO,Q>.iv)] = s^u i Ki -v)= s^JiC/v)- 

Turning now to the nonlocal component, (and neglecting the phase factor a'^ in q^), 
we have 

A^ «!y,;y + l A^ 

< loga^,^+ia^+i,^ - loga^^a^+i,^+i,loga,,+i,,, - logo,,,,,+i > . 
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From (2.14) 

< loga^,^+ia^+i,^ - loga^^a^+i,^+i,logai,+i,^ > (3.2) 

(-1)- 



n'^{^r]) 



n-l 



Qfj.,u+i ~ Qfi+i,u+i — Q^,,v + Q^,,u — Q^,+i,v+i + Q^,+i,u 







and the same holds for the quantity < ■ ,logajy,,y+i >. Hence < P^iiO^Qj^iv) >= 
and the proof of Theorem 3.1 is complete. 

The last computation shows that for all n, u, 

<P^i,ai,^^j^i>=<p^,a^j^i^^>=Q. (3.3) 

This fact will be used in §5. 

4. Action-angle variables; self-adjoint case. 

The canonical coordinates constructed in §3 are complex; but when L is Hermit- 
ian symmetric, the scattering data possesses an additional symmetry which allows 
us to construct real canonical coordinates. In the following, R denotes the order 
reversing matrix 

n 

R = 2_^ ^k,n-k+l 
k=l 

corresponding to the permutation (1, n)(2, n — 1) . . . ; as before Cjk is the matrix 
with a 1 in the jk place and zeroes elsewhere. 

Lemma 4.1 There is a strictly lower triangular matrix G(x), G(x) ^0 as \x\ -^ 
oo, such that 

ruL* (x, z) = {I + e{x))nz''-^R[mL{x, z)-^]*RJ{zy (4.1) 

for any z G C\E such that m{x, z) exists. 

In this lemma the two wave functions are given in their local representations. 
When n is even, z = a^z for some k; while if n is odd, z = —a^z for some z. 
Therefore, by the invariance of the wave functions under the action z ^ az of 
Lemma 1.2 we may write (4.1) as 

rriL^x, z) = {I + e{x))n z''-^R[mLix, a"" (z))]* RJ (z)* (4.2) 

where 

a(z) = —z. 
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This expresses the relationship between the wave functions for L and L* at points in 
the fundamental domain. In [BDT] this identity is used to establish symmetries of 
the jump data when L is self-adjoint. We modify this result here to make it directly 
applicable to the form of the scattering data (1.12). We prove the following: 

Theorem 4.2. Let ^ belong to the ray T^j and $, to T^k- Then 

RaL^OR = 7v^J-'[aL{0~'rJj^j- (4.3) 

The matrices aL and ql* in (4.3) are in their local representations; ttj is the 
permutation matrix associated with the ray Ej; and Jj = diag(aj, . . . ctn) in the j^^ 
ordering. We first need to prove the analogue of Lemma 4.1 for the wave functions 
m{x, ^) defined for ^ in E. We indicate the necessary modification of the proof of 
Theorem 9.5 in [BDT, pp. 41, 42]. When ^ belongs to Ej, ^^(x,^) is the unique 
solution of 

Dm = J^m — ^mJj + qm, 

lim m = (i(^)Aj, sup |m(a;,^)| < +oo. 

X — » — OO J. 

Lemma 4.3. Let ^ G aj, and ^ G E^. Then 

(i) TTfcJfcTTfc = Jfc + i; Afc + lTTfe = Afc 

(ii) RJjR = Jjt+i; RAjR = Ak+iJk+i = AkJki^k 
(iii) RtTjR = TTk] RU^R = %. 

Proof: The first statement in (i) is a consequence of the defining relation of the 
permutations ttj and of the ordering of the roots of unity; the second statement is 
a result of the fact that only the columns of A are permuted in changing represen- 
tations. From (1.6), 

Re idnZ > Re idn-iz > ■ ■ ■ > Re idiz for ^ G O 

and it follows that complex conjugation reverses the ordering of the roots. Hence 

RJ{z)R = J{z)\ RU^R = n^, and Rn^R = tt^. 

This establishes (iii). A direct calculation shows that RA{z)R = A{z)J{z). As 
z ^ ^ G Ej from Oj,z ^ ^ G Ejt from O^+i; so the identities in (ii) follow. This 
completes the proof of Lemma 4.3. 

Repeating the derivation of (9.14, [BDT]) with the wave function 
m{xj$,) we find that 

g{x) = nC-'R[mL{x,0-'rRJk+i 
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satisfies 

Dg = J^g - ^gJk+i + Rq*Rg, g{x) -^ A^+i as a; ^ -cx). 

Hence w = givk satisfies 

Dw = J^w — ^wJk + Rq*Rw, w ^ Afc as x ^ — oo. 

As in Tfieorem 9.5 of [BDT], tliere exists a matrix function G{x) = (/ + 0(a;)) sucli 
that 

mL*{x.O = Gix)nC~'R[mL{x,0~'rRJ*k+i^k. (4.5) 

Proof of Theorem 4.2: From (1.12) we have 

= hm H.-e-^^^'^'=Ari$^.(x,0 

X — >CXD ^ ? 

where 

$L* (x, = niL* (x: Oe'""^'" and A^- = di^Ak- 

From (4.5), 

aL*{0= lim U,e'''^^^Aj^G{x)nC~^R[mL{x,0~^]*RJl,^nke'''^^^. 

X — >oo ^ ? 

= R hm U,njJ*{e'''^'^^yA,*[m{x,0~^]*Jj{e-'''^'^^y7ijR 

X — »oo -' 

by repeated use of Lemma 4.3. Since tcj and H^ have the same block structure, they 
commute, so 



j^j- 



-' X — >CXD 

Now \l/$ = a and \1/ = ^A7^ so a~^il) = e~'^^^'^^ m^^ A^; therefore 

hm U.e-'"'^-^^ ip ^ 1 

X — »oo 

and Theorem 4.2 foUows. 

Complex conjugation does not preserve EqUEi, where the canonical coordinates 
were defined in §3, but we may combine it with the invariance under rotations of 
Lemma 2.1 to obtain an involution. For ^ in Eq U Ei, let ^* be the unique point 
in Eq U El which is on the same orbit as ^ under rotations by a. When n is even. 
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^* = ^; while when n is odd, the involution interchanges Eq and Ei. In both cases 

in = 1^1- 



Theorem 4.4. For L = L* : 



1 -^^-r 1 



«ii(0 = 77^, A.(0 



(4.6) 



Proof: The first identity follows directly from (4.3). To prove the remaining iden- 
tities we use the relation 



Bn-AC) ^ 



A.(0 






(6, (4.7) 



which is a block by block restatement of (4.3). In deriving (4.7) recall that Jq = 
diag(l, a, a~^, . . . ) and Ji = diag(a, 1, a^, a~^, . . . ); we work in the O*'^ representa- 
tion. 

The second identity in (4.6) is immediate from (4.7). We leave it to the reader 
to verify the third; while for that involving the Qj^'s, we have 



Qn-ui^ = 1— log '■ a" '^ = — -loga''^ a 

Am an-u,n-u+i 47rz a "a^+i 



n — v 

'Jt 

V 



^-—\o^o"— = -— -loga ""- '- = qAi). 

Am a^+i^^ Am a^^^+i 

These symmetries allow us to construct real canonical variables in the self- 
adjoint case. For odd n they are 

V2RepA^), V2ReqA0. 

V2ImpA0. V2ImqAC): 1 < ly < ■ ■ ■ < n/2. 

When n = 2m is even, pm and g^ are real and these variables are to be added to 
the above list. These real variables also satisfy the canonical commutation relations 
when ^ and r] lie on the rays Eq and Ei. For example when fj., u < n/2, 

{V-2 Re pAO. ^2 Re ,,(„)} = { ^-«' +^-''^'' . *iM±|zi^} 
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since d^^n-n = when /U, z/ < n/2. 

5. Complete Integrability of the Gel'fand-Dikii Flows. 

We may now demonstrate the complete integrabihty of the Gel'fand-Dikii flows. 
Lax pairs for the flows (1.1) can be obtained [BDT] as a pair of flrst order operators 
Dj: = D — zJ — Qz and Dk,t = Dt — (zJ)'' — Bn,k, where Dt = —id/dt, and the 
i?n,fc are n x n matrices whose entries are polynomials in z^uq, . . .Un-i and their 
derivatives up to order k — 1. Under the flows given by the zero curvature condition 
[Dx, Dk,t] = 0, the scattering matrix a evolves according to 

a(e,t) = e^(«^)'*a(e,0)e-^(«'^)'*. (5.1) 

As in the case of integrable systems built on flrst order n x n isospectral operators, 
the scattering transform decouples the modes, in the sense that the evolution of 
the scattering matrix at ^ G E is independent of ^' 7^ ^. It is easily seen that the 
pointwise flows (5.1) are Hamiltonian with respect to the pointwise Poisson bracket 

Til 

(ajk, aim) = —-——aikajm[sgn{l - j) + sgn{k - m)]. (5.2) 

n{—c,) 

Relative to the bracket (5.2), the local component of (2.14) is (ajfc(^), aim{'fl))^{C/v)j 
since ^ 6{^ - i]) =d{^/r]). 

In what sense do the flows (5.1) comprise a "complete set of commuting Hamil- 
tonian flows?" For flxed ^ G E the minimal scattering data at ^ is given by the 
upper and lower triangular factors v'^. of B,^ with diag v^ = 1 (cf. [EC] or [BSl]). 
Each v^ has one nontrivial entry, so the scattering data on E consists of a total 
of 2(n — 1) functions deflned on the two half lines Eq and Ei. (Compare with the 
n — 1 functions -uq, . . ■Un-2 deflned on the real line). There are [n/2] degrees of 
freedom for each ^ G Ei and n — [n/2] — 1 degrees of freedom at each point ^ G Eq. 
The question of complete integrability of the Gel'fand-Dikii flows is now reduced to 
a flnite dimensional problem, namely that of flnding n — 1 pointwise independent 
commuting flows on the scattering data for each flxed ^ G E. It is clear that, for 
flxed ^ G E, the flows (5.1) comprise just such a maximal set of commuting flows. 

The Hamiltonians for the flows (5.1) are generated by a linear combination of 
the Piy's. Straightforward computations using the properties of the brackets (e.g. 
the derivation property) show 

{Pu,a^^^+i) = -2TciSiyf^a^^^+i, 

{p^, a^+i,^) = 27rz5i.^a^+i,^, {p^y, a^^) = 0. 
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Therefore the flow deflned by 

is given by Biy{t) = ex.p{iaiyt}Biy{0) ex.p{—ia^t} and B^{t) = -B^(O) if /U 7^ z/, where 
Sjy is the 2x2 diagonal block and a^, is the diagonal matrix with a 1/2 in the z/*'^ 
entry and -1/2 in the (z/ + 1)** entry. 

For ^ in Ej, let av± = a± (assuming this factorization exists), where v+, a_ 
{v-, a_|_) are lower (upper) triangular in the j*'^ representation; and diag v± = 1. 
The matrix v{^) introduced in §1 as the scattering data equals vZ V-\- ([BC]). Define 
d± =diag a±. By direct computation (cf. Lemma 4.2 of [BSl]) we obtain 



d_^6+ = diagf — 
Van 



Ai aiia22 A2 



' A ' ■ 

0^22 ^1 «33ti44 



on the odd numbered rays, with a corresponding result holding on the even num- 
bered rays. By (3.1) 

log 51^5+ = — -diag(pi,-pi,p3,-p3,...) (5.3) 



27ri 2-1X1 

on the odd rays, again with a corresponding result on the even rays. Therefore the 
pointwise fiows (5.1) are generated by 



2iii 
The fiows (5.1) are generated by the direct integral of these pointwise Hamiltonians 



^k = ^ [ (-0"tr[(eJ)^ log 51^5+] f 

2m Jj2 C 



(5.4) 



Indeed, consider the functionals 



Fi^ = J aim{^,t)f{0- 

and the associated fiows defined by DtFim = {Fim^'Hk}- By (3.3) the nonlocal part 
of these brackets vanishes, and by (5.2) and (5.3), 
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Since f is an arbitrary element of T>{'E), 

and the corresponding flows are precisely the Gel'fand-Dikii flows (5.1). 

These integrals are well defined if the potential q is in the Schwartz class; for 
then ([BDT], Corollary 5.18) log5Z^5_|_ tends to zero rapidly as ^ tends to infinity 
along the rays of E. Moreover ([BC]) the diagonal matrices 6± are the boundary 
values on E of the sectionally meromorphic function d{z) obtained as 

5{z) = lim (j){x,z)e-''''^^'\ 

X — >oo 

In the absence of discrete components of the scattering data, 5{z) is holomorphic 
in 0± and its entries have no zeroes. 

For integrable systems based on n x n first order isospectral operators an 
(n — l)-parameter family of generating functions for the Hamiltonians is given by 
tr[//log5(2)], n a traceless diagonal matrix ([S],[BC]). We have a similar result here, 
with a single generator. 

Theorem 5.1 The generating function for the Gel'fand-Dikii flows is 

n{z) = {-iy+h''tr[{I + J + J2 + ■ ■ ■ + J^-i) log6{z)]. 

Note that n~^{I + J + ■ ■ ■ + J"~^) projects onto the null space of J — I. The 
generating function has branch points at the zeroes of dj (z) , which correspond to the 
discrete part of the scattering data; but we again omit this case and prove Theorem 
5.1 under the assumption that the scattering data has no such discrete part. Under 
those assumptions H(^) is single valued and analytic in each component of C\E. 

Proof: The sectionally analytic function ^k{z) = tr[j'^(z)logd{z)] is 0{z~^) as 
2 ^ oo. Its jump across E is [^k] = tr[J(^) log6Z^S-\-] (using the same representa- 
tion for 6- and 6+ on E) so by the Plemelj formulas 

If tr[j'-(01ogC'a+] ,^ 

*'('-' = 2^7, — r^z — <^- 

(There is a weak singularity at ^ = 0: d{z) behaves like a power of 2; as 2; tends to 
zero in the closure of any of the sectors 0,j [BDT] ; so log 5Z 5_(_ has a logarithmic 
singularity at the origin.) Since zJ{z) is invariant under z — > az, it follows that 
$fc(a2) = a~''^k{z): and 



zm Jy, "~ t, — d^ z 
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From the identity 






n- 



which may be verified by comparing the residues of both sides with respect to the 
variable ^, it foUows that 



J=0 

oo 






oo ^ 

^ ' / J ynj-\-k 



i=-i 
The result is obtained by summing from k = to n — 1. 

We would like to thank Boris Konopelchenko for helpful discussions during the 
preparation of this paper. 
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